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Abstract 


This paper reports a nonlinear fuzzy modeling study of a molten carbonate fuel cell (MCFC) stack by an identification method. MCFC is a 
complex nonlinear, multi-input and multi-output (MIMO) system that is hard to model by traditional methodologies. The Takagi-Sugeno (T-S) 
fuzzy model is suitable to model a large class of nonlinear MIMO system. In this paper, a MIMO T-S fuzzy model is used to represent MCFC. An 
identification method is used to determine both the nonlinear parameters of the antecedents and the linear parameters of the rules consequent in the 
T-S fuzzy model. The simulation tests reveal that obtained T-S fuzzy model using the identification method can efficiently approximate the static 
and dynamic behavior of a MCFC stack. Furthermore, based on this proposed T-S fuzzy model, valid control strategy studies such as predictive 


control, robust control can be developed. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


The molten carbonate fuel cell (MCFC) is an energy conver- 
sion device that produces electricity by the indirect combination 
of hydrogen and oxygen to water via an electron-carrying elec- 
trolyte. MCFC is one of the fuel cell technologies that have 
proven efficiency and environmental performance [1]. In the 
development of fuel cell technologies, the mathematical model- 
ing is an important tool, which has the capability of predicting 
the fuel cell performance [2]. The results obtained from a reliable 
and effective model can be very useful to guide future research 
for fuel cell improvements and optimization. At the same time an 
effective model of the MCFC stack, is a prerequisite for control 
analysis and controller design. 

As far as we known, the MCFC is a nonlinear multi-input 
and multi-output system. It is very difficult to model the MCFC 
system by the traditional methods. During the last several 
decades, various mathematical models have been established 
in the research on the internal mechanisms, ranging from a one- 
dimensional model to a three dimensional model [3-5]. These 
models are very useful for cell design and performance analy- 
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sis. However, these models are too complicated to be used for 
control studies. 

To meet the demands of developing valid control strategies, 
some researchers have attempted to establish novel fuel cell 
models by statistical data-driven approach. The artificial neu- 
ral network (ANN) has been used to derive a solid oxide fuel 
cell (SOFC) model from the experimental data quickly [2]. Shen 
et al. [6] proposed a RBF neural networks identification technol- 
ogy to set up the nonlinear temperature model of MCFC stack. 
Jemei et al. [7] utilized ANN methodology for proton exchange 
membrane fuel cell (PEMFC) model. The AutoRegression with 
exogenous signal (ARX) identification algorithm was applied 
to computer linear SOFC system models [8]. The Hammerstein 
nonlinear system approach was used for identification of the 
SOFC model [9]. The support vector machine (SVM) was used 
to modeling study of the PEMFC and SOFC system [10,11]. 
However, most of modeling methods are in the point of steady- 
state models (—V characteristics) or the single output, and these 
black-box models do not give much insight into the system being 
modeled. 

Fuzzy models have proved to be useful in nonlinear dynamic 
system modeling. Several modeling methods based on fuzzy 
reasoning have been proposed in recent years, which can be 
roughly classified into the fuzzy relational model, the neural- 
network-based fuzzy model, the T-S fuzzy model and the 
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Nomenclature 


AR frequent factor of the resistance 

AR,, Ap, Ag coefficients of the cell inherent perfor- 
mance 

Cac) total concentration of anode (cathode) 

C i heat capacity of the cell unit (J (kg K)~!) 


E thermodynamic voltage (V) 
E? reversible potential of the cell (V) 
F Faraday constant (96,485 C mol!) 


Pars) anode (cathode) total inlet mole flow (mol s~!) 

PAS) molar flow rates of the ith reactant at the cell 
input of anode (cathode) 

Fit) molar flow rates of the ith reactant at the cell 
output of anode (cathode) 

F, (Fe ;) reaction rate of the ith reactant of anode (cath- 
ode) 

hac) channel height of anode (cathode) (m) 

hin hi) anode (cathode) inlet partial molar enthalpies 
(Jmol!) 

h? partial molar enthalpies at stack temperature 
(Jmol!) 

AHr, AH, activation energy of internal resistance, and 
AHo,, AHco, anode and cathode gas 


I current (A) 

J current density (A m~?) 

Ms mass of the cell unit (kg) 

Neell cell number 

p gas pressure 

Pac stack dc power (W) 

R gas constant (8.3145 J (mol K~!)) 


Ra i(Rc,i) anode (cathode) total rate of production of 
species (mol s7!) 


Rohmic ohmic resistance (Q m~?) 
sS cell active area (m?) 
T stack solid average temperature (exit temperature) 


(K) 

rT) inlet temperature of anode (cathode) (K) 

Uf(Uox) fuel (oxidant) utilization 

Vae cell voltage (V) 

ae) anode (cathode) inlet mole fractions 

X,;(%e;) anode (cathode) outlet mole fractions 

Za(Zc) impedance for electrode polarization of anode 
(cathode) (Q m7?) 


fuzzy basis function based model. Compared with other mod- 
els, the T-S model needs less rules, each rule’s consequence 
with linear function can describe the input-output mapping in 
a large range, and the fuzzy implication used in the model 
is also simple [12]. Takagi-Sugeno (T-S) fuzzy model has 
the ability to approximate a large class of static and dynam- 
ical MIMO nonlinear systems. There exist a large number 
of research papers on the application of model identifica- 
tion successfully and control strategies for T-S fuzzy model 
[13-15]. However, the concrete study of modeling multivari- 


ate MCFC using T-S fuzzy model cannot be found in prior 
papers. 

In this paper, a nonlinear MIMO T-S fuzzy modeling of a 
molten carbonate fuel cell (MCFC) stack is built with an iden- 
tification method. T-S fuzzy model consists of if-then rules 
with fuzzy antecedents and mathematical functions in the con- 
sequent part. Antecedent identification by which the MCFC 
inputs—outputs variable space is divided into many subspaces 
(rules) is implemented based on the principle of Fuzzy C-Means 
clustering method. In every rule, the consequent part parame- 
ters that make the nonlinear MCFC characteristic to be fitted 
by a linear model are identified by using the Kalman filtering 
algorithm. In the process of fuzzy modeling, a novel way is intro- 
duced to determine satisfactory number of rules. In our study, the 
proposed MIMO T-S fuzzy model of MCFC is built with data 
obtained from a physical model of 10 kW MCFC. The encourag- 
ing results from this study demonstrate that there is a potential to 
introduce the nonlinear fuzzy modeling to development of fuel 
cells by the data obtained from various fuel cell experiments in 
the future application. 


2. MCFC dynamic physical model 
The proposed model is based on the following assumptions: 


(1) Stream mixture thermodynamic properties follow the ideal 
gas mixture law. 

(2) A uniform gas distribution among cells is supposed. 

(3) All cells have the same temperature and current density. 

(4) Fuel processor dynamics are not considered. 


2.1. Electrochemical model 


The basic structure of a single cell is shown in Fig. 1. There 
are two main sets of chemical kinetics. This chemical kinetics 
are those associated with anode reactions, cathode reactions. 
In the following equations, R vector is used to represent rate 
(respectively) of five individual species, which are hydrogen, 
carbon dioxide, water, oxygen and nitrogen. 


e anode: Hz + CO32~ > CO? +H20 + 2e7 


+ 


R =| Neen! Neend Neen! 
j 2F 2F 2F 


Current 
collector 


Electric load 


1/20,+CO, OQ= 


Fig. 1. Principle of electricity generating in a MCFC single cell. 
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where using the ideal gas law, the exit stream total concentration 
becomes Ca(c) = Pa(c)/RT°. The compartment volume Vac) equals 
to sha). The fuel utilization up = Nee Js/2 FF; in y x . 


2.2. Thermal model 


The fuel cell power output is closely related to the temper- 
ature of the cell unit. Assuming that energy accumulates only 
in the stack solid mass, gas mixtures are ideal, and exit stream 
temperatures are equal to the solid stack temperature, the energy 
conservation equation is [16]: 


5 
_ pe a q,- oi - YOR Rai 
i=1 


va > 


+ F} 


Ohoss = Pac (4) 


5 
aD Niy- 
i=l 


Under the ideal gas supposition, the partial molar enthalpies are 
calculated by the formulation: 


T 
hi = href +f Cp,i(u) du (5) 
Tref 


And the coefficients of the specific heats cp, are encountered in 
standard reference tables: 


Cpi = ai + biT + ciT? + diT? + eiT" (6) 
2.3. Operating cell voltage 


The empirical relationships developed in previously reported 
work [17] are employed to estimate a voltage drop in a fuel cell. 


The cell-unit performance of the electric power generation is 
expressed: 


Vac = E — (Rohmic — Za — Ze) J (7) 
where E is given by the Nernst potential: 


1/2 
RT In Ptba Pd, PCO», 


E=E°+ (8) 
2F PH20,aPCO>,a 
and 
o _ 4184 x [58.3 — (0.0113 + 9.6 x 10-7T)T] 
B= (9) 
2F 
Raced Anr (10) 
ohmic = AR EXP | — RT 
AM,\ o5 
Za = Ar,T exp (- RT ) Pr. (11) 
AHo, \ -0.75 0.5 
Ze = ApT exp (ee) a Po, PCO 
AHco, \ -1.0 
+ AgT exp (- RT ) Peo (12) 


In Eqs. (8)-(12), T is the arithmetic average of cathode 
inlet and exit temperature. Anode and cathode partial pres- 
sures are arithmetic averages of inlet and exit gas partial 
pressures: 

Tin + T° 
T= ® 13 
5 (13) 
Xalc),i + Xa(o),i 


(14) 
2 Pa(c)/ Patm 


Pa(c),i = 


3. T-S fuzzy modeling 
3.1. Problem formulation 


We consider a MCFC system G(U, Y) as a MIMO sys- 
tem, let u(u € U C R’) be input variable and y(ye Y C R3) 
output variable. For the MIMO system, it can be divided into 
q multi-input and single output (MISO) system, and each of 
MISO can be fit by a fuzzy T-S model independently [14]. We 
define: 


(uO = ujk), ujk- mu, D] G=1,...,, 
Diko” = Dik), ..., wk — nu +D] G=1,....9 


with ny, j and ny, are the order of uj and y;, respectively, then 
each MISO subsystem can be denoted as follows: 


wk + 1)= fk) @C=1,...,9 (5) 


with x(k) = [ur Oo -a rD > Lyi D}o", - 
{yg(K)}g7] = = [Xki>---;Xk,] is the regression data vec- 
tor consisting of input-output data at the kth instant and 
before. 
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The T-S fuzzy model with linear consequents employed to 
fit the MISO subsystem in this paper is a collection of fuzzy 
rules, which is in the form of “If. . . then . . .” The ith rule of the 
lth output 37, ;(k + 1) is given by: 

Rii: If x(k) is Al, then 


Silk + D= Phot Dix, to + Pik G=L....0 
(16) 


where c is the number of rules, and A! = {Al TE 


Pr an is 
the set of membership functions associated to the ith rule. p! = 
[ Plo» Pas foe Phyl is the parameter vector of the ith submodel 
(rule). 

If the method of product fuzzy inference and weighted mean 
is employed, the output of the T-S fuzzy model is inferred as 


follows: 


Dein Milk) ilk + 1) 


k+ l= À 
ĵi ) SE u 


a7) 


milk) = | [u (18) 
j=l 


where uij(k) is the match degree of the component xg; with 
respect to At p- Lik) is the match degree of x(k) with respect 
to Al, namely the membership degree of x(k) in the ith 
tule. 

The antecedent part of “If...” in the form of fuzzy set, is 
equal to fuzzy partition of data space of the MCFC plant, and the 
consequent part of “then...” represented by a linear functional 
relation, is a linear composition of input—output variables. The 
identification problem that we want to solve is to find the fuzzy 
MCEC model of reasonable complexity that minimizes the error 
variance between the fuzzy MCFC model output and the real 
MCEC stack output. In this model, what is to be partitioned is the 
whole space spanned by regression data vectors. We separate the 
identification of the fuzzy model into two parts—antecedent and 
consequent identification. All the subsystems (y;(k + 1) = f(x(k)), 
l=1,...,q) can be identified independently in the same 
way, so the subscription / is omitted in the following 
description. 


{xk: k=1,...,N} be a set of N sample data from the real plant. 
Each sample xk = [xx,, Xk), - - - , Xk, | has n components. 

The FCM algorithm was introduced by Bezdek [18]. In this 
paper, FCM algorithm is used to partition the collection of N 
data from real MCFC stack into c fuzzy clusters. In the fuzzy 
modeling, the c means the number of fuzzy rules. The idea of 
FCM is using the weights that minimize the total weighted mean- 
square error: 


N c N c¢ 
Ju, V) = Y u" dk = SO ui) lle — vill? 


k=1i=1 k=1i=1 
(19) 


where mik is the membership degree of the kth sample in 
the ith cluster. u is the matrix «~={;%} which satisfies: 
Lik € [0, 1], Sk = |. The exponent m is a real number 
greater than unity. V = [v1, v2,..., Uc] represents the center of 
each cluster. The djz is the distance from xz to v;. Local minimiza- 
tion of the function J is accomplished by repeatedly adjusting 
the values of uix and v; according to the following relation: 


1 
ja Giz /d jy 


N m 
E Žr- Mi k*k 
po N m 
ee ik 


(20) 


Mik = 


(21) 


3.3. Consequent then-part identification 


The consequent part of the fuzzy rule is identified by using 
the Kalman filtering algorithm. 

Given the data x, where k represents the kth sampling, the 
(k+1)th predictive output $(k+ 1) of fuzzy model, can be 
obtained by formulation (17). Let By = ti(k)/ De pj(k), then 
(k + 1) can be written as: 


c 
Ik + D) = XO Bri (Pio + PiX H+ Pip Xky) (22) 
i 


While xk = [Xk,,Xk,---, Xk, ] (k=1, 2,...,N) represent the 
data of n dimension input variables of MCFC, and Y=[y1, 
y2, - . yy] are the outputs, and the k = (Bx,, Bko» -< , Be.) can 
be calculated. The outputs can be expressed from (22) using the 
matrix form as follows: 


Yon (23) 
Bir... Bree». XB... xu Ble-..  Xinßi1---  Xinßic 
Bar... Bre x21Po1-.. Xabre = Xan B21... Xa Bre 
h= (24) 
Bni--- Bne--» XniBni--- XNn1BNe-.- XNnBN1--- XNnBNed yxen+1) 
P = [P10, P20, - -< , PcOs P11, P21, +++» Pcl + +++ Pins 
3.2. Antecedent if-part identification Dimens Beal ii (25) 


Antecedent identification is implemented by fuzzy clustering 
based on the principle of Fuzzy C-Means (FCM) algorithm. Let 


The consequent parameter vector P can be estimate by Eq. 
(23). To void the computationally expensive matrix inversion 
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operation, the Kalman filtering algorithm is applied in this 
paper. The Kalman filtering algorithm enables calculation of 
the new, adapted values of parameters of the antecedent part, 
on the basis of the new sample and the known parameter values 
[19]. Here, we apply it to calculate the parameter vector P as 
follows: 


Skt thig Ong — hrti Pe) 


Peyi = Pet (26) 
ss Q + hes iSkhyy 
Shy heer, 
Sit = Sp— ee (27) 
Q + hk+1Skħk41 
k=0,1,2,...,N—1 (28) 


where Py is the estimated parameter vector. Po is zero vector. The 
hg+1 is the (k+ 1)th row of matrix h, and yg+1 is the real output 
of (k+ 1)th sample point. Sọ = a/ (Iis an identity matrix and a is 
a large positive number), and Q=exp(—WN/z) (N represents the 
iteration counter and z is a constant). 


3.4. Fuzzy modeling procedure 


The number of fuzzy rules is an important factor that affects 
the performance of a fuzzy model. The choice of an appropri- 
ate number of rules is crucial to the design of fuzzy systems 
because practical considerations usually need to keep the num- 
ber of rules as low as possible to reduce complexity. While 
too many redundant rules result in a complex fuzzy model and 
increase difficulties of implement, too few rules produce the 
fuzzy model cannot approximate a real plant. In the paper, the 
cluster number c in FCM algorithm corresponds to the number of 
“If... then” rules. However, for many kinds of clustering meth- 
ods, including FCM algorithm, clustering number c is always 
needed in advance. Here, a novel way [15] to determine the 
cluster number (namely the number of fuzzy rules) is intro- 
duced to the fuzzy modeling. Let the clustering method start 
with c=2. Then determine whether a new cluster center should 
be increased or not by evaluating modeling results. Root mean 
squared error (RMSE) is employed here to evaluate modeling 
results: 


RMSE = (29) 


where N is the number of sample data from the real MCFC stack, 
$x is the predictive output of fuzzy model and yx is the output of 
the real MCFC stack. If the modeling result is not satisfied, from 
the given sample data, a sample x; that is most different from 
the existing cluster centers V = [v1, v2, ..., Vc] can be fined as 
a new center v-+1. Where the following definition for ve+1: 


Ver = xe | k= argemin XO (u-n) (30) 


l<ij<c 


i$j 


Initialize partition matrix uo (randomly), 
fuzzyfication exponent m, and the stop 
criterion €(€ > 0) 


Estimate the Antecedent model 
parameters using FCM, by matrix uo 


Identify the consequence model 
parameters using Kalman filter algorithm, 


find outa sample Xk 
k=arg,min >, (Kik Mj) 
ISi j<c 
iej 


Update po = Litile+ien Using (20) 


Va = Q V =[vp Vy aos 


Fig. 2. Flowchart of the fuzzy modeling method. 


Start with V = [v1, v2,..., Vc+1] as initial cluster centers, 
and compute the new not-random partition matrix jo. Then 
FCM algorithm is applied to divide the sample data into c+ 1 
parts again. Do the above steps again until the result is sat- 
isfactory. Then the T-S fuzzy model can be obtained by the 
antecedent and consequent identification methods. The basic 
design steps of the modeling method are depicted in Fig. 2. 


4. Results 


As we know, output voltage and the temperature are the core 
of any fuel cell modeling. For a given MIMO MCFC stack, the 
output voltage and the temperature are influenced by many oper- 
ating parameters such as current density, temperature, pressure, 
fuel utilization (uf), oxidant utilization (uox), and gas flows, etc. 
However, in practice, the gas flows must vary within the allow- 
able gas utilization range considering the safe operating areas 
of a plant. Up to now, no model has ever been able to accom- 
modate all these operating parameters. Our T-S model is no 
exception. In our experiment, the MCFC can be regarded as a 
system with two-input (uf, uox) and two-output (Vac, T°), and 
current density J is considered to be a disturbance and other 
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operating parameters are held constant. Define the following 
input—output relations: 


Ivalk + 1) = fuk), Uoxlk), J(k), Yva (k)) (31) 
Sro(k + 1) = fulk), uox(k), J(k), yro(k)) (32) 


where Jy,., jo are the predictive outputs of fuzzy model and 
YVae> Yre are the outputs from the real MCFC stack. The order 
of inputs and outputs can be increased according to a concrete 
system. This above simplification does not impair the valid- 
ity of our study, which is aimed at modeling MCFC by the 
proposed method in this paper. The proposed fuzzy modeling 
method is used to identify a T-S fuzzy model of MCFC based 
on simulation data. 


4.1. Preparing simulation data 


The 10kW MCFC stack is used in the simulation. The stack 
consists of 25 cells with the anode and cathode in co-flow. The 
effective electrolyte area of the cell is 0.4 m7, and the electrolyte 
tile between anode and cathode consists of mixed carbonate of 
Li2CO3 and K2CO3. The compositions of fuel and oxidant are 
set at constant values. The dynamic physical model as shown in 
Section 2 replaces the real MCFC stack to generate the simula- 
tion data required for the identification of the T-S fuzzy model. 
The parameters of this fuel cell are given in Table 1. The data 
sources blocks developed in MATLAB based on Section 2 is 
shown in Fig. 3. For the purpose of identification, the dynamic 
physical model is excited with uniformly random input signals 
included the fuel utilization (40-90%), the oxidant utilization 
(10-80%), and the current density (100-1800 A m2). To obtain 
values at integer time points, the fourth-order Runge-Kutta 
method was used to find the numerical solution to the dynamic 
physical model in the simulation. A set of 10,000 data was col- 
lected from the simulation. The first 6000 data were used for the 


Table 1 

Parameters of the MCFC stack used in the fuzzy modeling 

Parameter Unit Value 

N cell 25 

Pac kW 10 

J Am”? 100-1800 
Vac Vv Variable 

ug 0.4-0.9 

Uox 0.1-0.8 

s m? 0.4 

ha m 1.2 x 107° 
he m 2x 10-3 
Cs kg m~? 7900 

AR 6.43 x 107? 
AHR kJ mol! —25.5 

AR, 8.11 x 107° 
AH, kJ mol`! _74.4 

Ap 2.10 x 10710 
AHo, kJ mol"! —83.4 

Ag 1.58 x 1075 
AHco, kJ mol"! —7.12 

qin K 873 

īm K 823 

x Mole fraction 0.64 

XCO Mole fraction 0.2 

ro Mole fraction 0.16 

Xoco Mole fraction 0.3 

XEN Mole fraction 0.553 

Xo Mole fraction 0.147 


identification of a T-S fuzzy model of the fuel cell, while the 
remaining 4000 data were used for validation purposes. 


4.2. Predicting with the MIMO T-S fuzzy model 


The parameters pre-specified in the modeling method as 
shown in Section 3.4 are the exponent m=2, stop criterion 


output voltage model 


Fig. 3. Data sources of T-S fuzzy modeling. 
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~ uf=70%,u0x=40%( physical model) 
` uf=40%,u0x=40%(physical model) 
21+ = uf=90%,uox=70%( physical model) 
me a uf=90%,uox=70%(T-S model) 
u ka uf=70%,uox=40%(T-S model) 
20 `u pA uf=40%,uox=40%(T-S model) 
`a E 
È 
~ 19F ES 
> “1 
E z “4 
ob ` 
S a. ee 
= 187 Ls pni 
© ` 
> x. 
a. 
17+ 1 
LS 
=- e 
" 
16 fF 7 
15 1 1 1 1 1 1 1 1 
0 200 400 600 800 1000 1200 1400 1600 1800 


Current density(A/m2) 


Fig. 4. J—Vac chart. 


£=0.15, and w=10?. After fuzzy modeling of Section 3.4, a 
T-S fuzzy model is obtained, which can be used to predict new 
input data. The estimated optimal number of fuzzy rules is equal 
to 12. The RMSE of output voltage obtained in training and test 
process is 0.1028 and 0.1172 respectively. The RMSE of temper- 
ature obtained in training and test process is 0.1249 and 0.1283 
respectively. The static J-Vaąc characteristics generated by the 
T-S fuzzy model showed good consistency with the physical 
model under various utilization, as can be seen in Fig. 4. Now 
the T-S fuzzy model is used to predict the dynamic characteris- 
tics of the physical model. The step changes in the stack current 
density (from 1000 A m~? to 1500 A m7?) and fuel utilization 
(from 80% to 85%) are applied with 50% oxidant utilization. The 
comparison between predicted (T-S model) and experimental 
(physical model) output voltage dynamic curves is represented 
in Fig. 5. At the same time, the predicted and experimental tem- 
peratures are shown in Fig. 6. From Figs. 4—6, we can see the 
obtained MIMO T-S fuzzy model can approximate the static 
and dynamic behavior of the physical MCFC model with good 
accuracy. 


— Physical model 
— T-S model 


Voltage(V) 


15 1 1 1 2 1 
0.7 0.8 0.9 1 1.1 1:2 3 1.4 1:5 1.6 


Time(s) 


Fig. 5. Output voltage dynamic response. 


930 
—— Physical model 

929 — T-S model 

928 

927 

926 

925 


924 


Temperature(K) 


923 


922 


921 


920 
0.7 0.8 0.9 1 1.1 12 1.3 1.4 LS 1.6 


Time(s) 


Fig. 6. Output temperature dynamic response. 


5. Conclusions 


To facilitate valid control strategy design and analysis of sys- 
tem stability, a fuzzy modeling study of MCFC is reported in this 
paper. An offline identification method for T-S fuzzy model is 
presented consists of the if-part identification and the then-part 
identification with a novel way to determine satisfactory num- 
ber of fuzzy rules. It is shown that the identified MIMO T-S 
fuzzy model of MCFC is more attractive in that it avoids using 
complicated differential equations to describe the stack, and the 
inputs—outputs static and dynamic characteristics of a MCFC 
stack can be predicted. In the future, based on this T-S fuzzy 
model, some control scheme studies such as predictive control 
and robust control can be developed. In addition, an online iden- 
tification algorithm for T-S fuzzy model with more variables can 
be considered. 
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